Control of spin coherence in semiconductor double quantum dots 
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We propose a scheme to manipulate the spin coherence in vertically coupled GaAs double quantum 
dots. Up to ten orders of magnitude variation of the spin relaxation and two orders of magnitude 
variation of the spin dephasing can be achieved by a small gate voltage applied vertically on the 
double dot. Specially, large variation of spin relaxation still exists at K. In the calculation, the 
equation-of-motion approach is applied to obtain the electron decoherence time and all the relevant 
spin decoherence mechanisms, such as the spin-orbit coupling together with the electron-bulk- 
phonon scattering, the direct spin-phonon coupling due to the phonon-induced strain, the hyperfine 
interaction and the second-order process of electron-phonon scattering combined with the hyperfine 
interaction, are included. The condition to obtain the large variations of spin coherence is also 
addressed. 

PACS numbers: 72.25.Rb, 73. 21. La, 71.70.Ej 
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I. INTRODUCTION 

The fast development of spintronics aims at making 
devices based on the electron spin. Semiconductor quan- 
tum dots (QDs) are one of the promising candidates for 
the implementation of quantum computations^ 2 -^ be- 
cause of the relative long spin coherence time, which has 
been proved both theoretically^ and experimentally! 7 ^^ 
Among different kinds of QDs, double quantum dot 
(DQD) system attracted much more attention recently as 
there is an additional coupling between two QDs in both 
verticali^iii^ and parallel^^^ DQDs, which can be 
controlled conveniently by a small gate voltage. There- 
fore spin devices based on DQDs can be designed with 
more flexibility. So far many elements of the spintronic 
device, such as quantum logical gateS ) 16 i 17 spin Alter a 18 i 19 
and spin pumps™ were proposed and/or demonstrated 
based on DQD system. Specially, in our previous work, 
a way to control spin relaxation time (T±) induced by the 
spin-orbit coupling (SOC) together with the electron- 
bulk-phonon (BP) scattering in DQDs by a small gate 
voltage was proposed^ However, according to our lat- 
est study^ the spin relaxation can be controlled by other 
mechanisms also, if calculated correctly. In the present 
paper, we include all the spin decoherence mechanisms 
following our latest study in the single QD system, and 
apply the equation-of-motion approach to study the spin 
decoherence in DQD system. By this approach, not only 
the spin relaxation time but also the spin dephasing time 
(T2) can be obtained. We show that both the spin relax- 
ation and the spin dephasing can be manipulated by a 
small gate voltage in DQD system. Especially, the large 
variation of spin relaxation still exists even at K. The 
DQD system studied here can be realized easily using the 
present available technology. 

We organize the paper as following. In Sec. II we set 
up the model and briefly introduce different spin deco- 
herence mechanisms. The equation-of-motion approach 



is also explained simply. Then in Sec. Ill we present our 
numerical results. We first show how the eigen energies 
and eigen wave functions vary with the bias field in Sec. 
IIIA. Then the electric field dependences of spin relax- 
ation and spin dephasing are shown in Sec. IIIB and IIIC, 
respectively. We conclude in Sec. IV. 



II. MODEL AND METHOD 

We consider a single electron spin in two vertically cou- 
pled QDs with a bias voltage Vd and an external magnetic 
field B applied along the growth direction (z-axis). Each 
QD is confined by a parabolic potential V c (r) = im*Wgr 2 
(therefore the effective dot diameter do = y/fvK /m*oJo) in 
the x-y plane in a quantum well of width d. The confining 
potential V z (z) along the ^-direction reads 



V z {z) 



eEz + \eVd, \a < \z : 

eEz+^eVd + Vo, \z\ < \a 
00, otherwise 



< ¥ 



(i) 



in which Vq represents the barrier height between the 
two dots, a stands for the inter-dot distance and E = 
Vd/(a + 2d) denotes the electric field due to the bias 
voltage Vd- Then the electron Hamiltonian reads H e = 
P 2 /(2m*) + V c (r) + V z (z) + H z + H SO , where m* is 
the electron effective mass and P = —ihV + |A with 
A = (B/2)(—y,x,0) being the vector potential. Hz = 
\g^sBa z is the Zeeman energy with g, and <x being 
the (/-factor of electron, Bohr magneton and Pauli ma- 
trix respectively. H so is the Hamiltonian of the SOC. In 
GaAs, when the quantum well width and the gate voltage 
along the growth direction are small, the Rashba SOC 21 
is unimportant Therefore, only the Dresselhaus term 2 ^ 
contributes to Hso- When the quantum well width is 
smaller than the QD diameter, the dominant term in the 
Dresselhaus SOC reads H so = i J2\ 7a( — Px&x + Py G y)-> 
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with 7^ = j Q (P 2 )/h 2 . 70 denotes the Dresselhaus co- 
efficient, A is the quantum number of z-direction and 
(P 2 ) x = -h 2 J r z ,x(z)d 2 /dz 2 ip z ^(z)dz, where V,,a (A = 
1, 2, 3 • ■ ■ ) is the eigen wave function of the electron along 
the z-direction^ The electron eigen energy and wave 
function in the cc-y-plane can be obtained by the exact 
diagonalization approach^ 

The interactions between the electron and the lattice 
lead to the electron spin decoherence. These interac- 
tions contain two parts, one is the hyperfine interaction 
between the electron and the nuclei, the other is the 
electron-phonon interaction which is further composed of 
the electron-BP interaction H ep , the direct spin-phonon 
coupling due to the phonon induced strain H str ain and 
the phonon-induced (/-factor fluctuation. We briefly sum- 
marize these spin decoherence mechanisms and the de- 
tailed expressions can be found in Ref. S (i) The SOC to- 
gether with the electron-BP scattering H ep . As the SOC 
mixes different spins, the electron-BP can induce spin re- 
laxation and spin dephasing. (ii) Direct spin-phonon cou- 
pling due to the phonon-induced stain H stra in ^ Because 
this mechanism mixes different spins and also is related to 
the electron-phonon interaction, it can induce spin deco- 
herence alone, (iii) The hyperfine interaction H e i*& It is 
noted that the hyperfine interaction alone only induces Ti 
since it only changes the electron spin, but not the elec- 
tron energy, (iv) The second-order process of hyperfine 
interaction combined with the electron-BP interaction 
Vef-ph = K2)[E ro ^ I ((^| J ffe P |m)(m|ff eJ (r)|4))/(^ 1 - 
em) + T. m ^.MUHei{v)\m){m\H ep \l 1 ))/{e t2 - e m )](h\, 
where \£i) (i = 1, 2, 3 • ■ • ) and are the eigen state and 
eigen energy of H e , respectively. Although the hyper- 
fine interaction can not induce T\ alone, it is noted that 
the second-order process, combined with the electron- 
BP interaction, can induce both T\ and Ti- The other 
mechanisms, including the first-order process of hyperfine 
interaction combined with the electron-BP scattering 27 
and the g-factor fluctuation 2 ^ have been proved to be 
negligible s 5 

Due to the SOC, all the states are impure spin states 
with different expectation values of the spin. For finite 
temperature, the electron is distributed over many states 
and therefore one has to average over all the involved 
processes to obtain the total spin relaxation time. The 
Fermi-Golden-Rule approach calculates the spin relax- 
ation from the initial state to the final one whose major- 
ity spin polarizations are opposite. However, the aver- 
age method is inadequate when many impure states are 
included. 5 Also the Fermi-Golden-Rule approach can not 
be used to calculate the spin dephasing time. Therefore, 
in this paper we adopt the equation-of-motion approach 
for many-level system with Born approximation devel- 
oped in Ref. [5j. When the spin dephasing induced by the 
hyperfine interaction is considered, as the slow relaxation 
of nuclear bath compared to the electron, the kinetics 
is non-Markovian. Moreover, because of the Born ap- 
proximation, this equation-of-motion approach can only 
be applied for strong magnetic field (> 3.5 T). 29 There- 



fore, the pair-correlation method^ is further adopted to 
calculate the hyperfine interaction induced T% for small 
magnetic field. What should be emphasized is that the 
SOC is always included, no matter which mechanism is 
considered. It has been shown that its effect to spin de- 
coherence can not be neglected. 5 
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FIG. 1: (Color on line) (a) The Lowest two eigen energies 
along the z-axis E\ (A = 1,2) and the energy difference 
A£i2 = £i — £i vs. the bias field E. Note the scale of A£i2 
is on the right side of the figure, (b) Square of the absolute 
value of the ground state wave function along the z-axis at 
two typical bias fields. In the calculation, the well width d = 5 
nm, the inter-dot distance a = 10 nm and the barrier height 
V = 0.4 eV. 



III. NUMERICAL RESULTS 

Following the method addressed above, we perform 
the numerical calculation in a typical vertically coupled 
GaAs DQD with barrier height Vb = 0.4 eV, inter-dot 
distance a — 10 nm and well width d = 5 nm. The 
GaAs material parameters and the parameters related to 
different mechanisms are the same with those in Ref. [H. 



A. Electric field dependence of eigen energy and 
eigen wave function along z-axis 

The eigen energy and eigen wave function along the 
z-direction are obtained numerically. In Fig. [T]Ja), the 
lowest two eigen energies £ i and £2 along the z-axis and 
their difference A£ 12 = £ 2 — £1 are plotted as functions 
of the electric field E. It can be seen that the energy dif- 
ference A£i2 increases quickly with the electric field E 
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when E is larger than 0.1 kV/cm. The eigen wave func- 
tion of the ground state along the z-axis also has large 
variation with E. It can be seen clearly in Fig. [ljb) that 
when E is very small (1.2 x 1CP 2 kV/cm), the wave func- 
tion of the ground state locates at the two wells almost 
equally. However, when E is large enough (89 kV/cm), 
the wave function of the ground state locates mostly at 
the quantum well with lower potential. The physics of 
such bias-voltage-induced quick change of eigen energy 
and eigen wave function can be understood as what fol- 
lows. Because of the large barrier height V and/or large 
inter-dot distance a, the two quantum dots are nearly in- 
dependent and the eigen wave function along the z-axis 
of the lowest subband spreads equally over the two QDs 
when the source-drain voltage is very small. Therefore at 
this time the energy difference between the lowest two en- 
ergy levels along the z-axis A£±2 is very small. However, 
with the increase of the source-drain voltage, electron 
can tunnel through the barrier and the wave function is 
almost located at one dot with lower potential and there- 
fore A£i2 increases quickly with E. 



B. Spin relaxation time Ti vs. Electric field E 

The spin relaxation due to various mechanisms at dif- 
ferent magnetic fields B and quantum dot diameters do 
are plotted as functions of electric field E in Fig. [2] 
The temperature T = 4 K. It is seen that with the in- 
crease of E, the spin relaxations induced by all the three 
mechanisms almost keep unchanged for small E (when 
E < 0.1 kV/cm), then increase a little and reach max- 
imums around 0.3 kV/cm. What is interesting is that 
when E is increased to around 1.1 kV/cm, the spin relax- 
ations are suppressed very quickly over a small window of 
E. Therefore, the total spin relaxation can be controlled 
effectively with a small value of the variation of the bias 
field AE. For example, in Fig. Ufa), the spin relaxation 
shows ten orders of magnitude variation when the elec- 
tric field E changes from 0.1 to 10 kV/cm. This can be 
understood as following. All the three mechanisms are 
related to the electron-BP scattering which is affected 
sensitively by the phonon wave length. The scattering 
becomes most efficient when the phonon wave length is 
comparable with the dot size. It is also known that the 
spin relaxations between the first and second subband 
are dominant for small electric fields Therefore when 
E < 0.1 kV/cm, as the energy difference between the 
lowest two subband along the z-direction A£i2 is almost 
constant (see Fig. [I] (a)), the spin relaxation keeps nearly 
unchanged. When the phonon wave length is comparable 
with the dot size, the electron-phonon scattering becomes 
most efficient. Therefore the spin relaxations show max- 
imums. However, with the further increase of energy dif- 
ference AS 12 by the bias voltage, the phonon wave length 
becomes larger than the dot size. Consequently the spin 
relaxation decreases very quickly over a small window of 
AE. 
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FIG. 2: (Color on line) T^ 1 induced by different mechanisms 
vs. the electric field at different magnetic fields and dot di- 
ameters. T = 4 K. Curves with ■: by SOC together with 
the electron-BP scattering; Curves with •: by the direct spin- 
phonon coupling due to the phonon-induced strain; Curves 
with T: by the second-order process of the hyperfine interac- 
tion together with the BP (V}p_ ph ). 



Now we focus on the variation magnitude of the spin 
relaxation with the bias field under different conditions. 
The largest variation of T\ (10 orders of magnitude) hap- 
pens at small magnetic field B = 0.1 T and small diam- 
eter dg = 10 nm. However, for larger d (20 nm in Fig. 
H£b)) or larger B (1 T in Fig. HJc)), the variations of the 
total spin relaxation decrease by several orders of mag- 
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nitude. It is further seen that when both do and B are 
increased [do = 20 nm and B = 1 T in Fig.[^d)], the vari- 
ations of the total spin relaxation is even smaller. This 
is because the spin relaxation induced by the electron- 

(3) 

BP interaction and h decreases with B and d in 

the high electric field region, where electron is mostly 
confined in one dot. This is similar to the single QD 
case£ Therefore to achieve large control of spin decoher- 
ence, the magnetic field B and dot diameter d should 
be small. 



10" 

7 a io- 2 



io- 4 



10- 




T = K 



io- 8L - 

IO -2 



10" 1 10° 
E (kV/cm) 



10' 



FIG. 3: (Color on line) 7\ 1 induced by different mechanisms 
vs. the electric field. In the calculation, do = 10 nm, B — 0.1 
T and T = K. Curve with ■: by SOC together with the 
electron-BP scattering; Curve with •: by the direct spin- 
phonon coupling due to the phonon-induced strain; Curve 
with ▼: by the second-order process of the hyperfine interac- 
tion together with the BP (v}p_ ph ). 



We further investigate the electric field dependence of 
spin relaxation at T — K and the results are shown in 
Fig. [3] It is interesting to see that the spin relaxation 
induced by the second-order process of hyperfine interac- 
tion combined with the clcctron-BP scattering (V^j_ ph ) 
still has large variation with E (5 orders of magnitude 
variation when E changes from 0.1 to 10 kV/cm). How- 
ever, the spin relaxations induced by the other two mech- 
anisms keep almost unchanged. This is because when 
T = 0K, the electron only locates at the lowest orbital 
level and the spin relaxation happens only between the 
lowest two Zeeman sublevels which keeps unchanged with 
E. Therefore the large variations of spin relaxations in- 
duced by the SOC together with the electron-BP scatter- 
ing and the strain-induced direct spin-phonon coupling 

(3) 

no longer exist. However, for V^j_ ph , which is the second- 
order process scattering, the middle states m) can be 
higher levels as the hyperfine interaction can couple the 
spin-opposite states in different subband along the z-axis. 
The energy differences between the middle states and the 
initial/final states increase with E and therefore the spin 

(3) 

relaxation induced by h decreases with E quickly. 
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FIG. 4: (Color on line) T-f 1 and T 2 _1 induced by different 
mechanisms vs. the electric field at different magnetic fields 
and dot diameters. T = 4 K. Curves with ■: by SOC to- 
gether with the electron-BP scattering; Curves with •: by 
the direct spin-phonon coupling due to the phonon-induced 
strain; Curves with T: by the second-order process of the 
hyperfine interaction together with the BP ); Curves 

with by the hyperfine interaction. 
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C. Spin dephasing time T2 vs. Electric field E 

Now we turn to study the variation of spin dephasing 
with the electric field E and the results under different 
conditions are summarized in Fig. [5] It is seen that the 
spin dephasing induced by the SOC together with the 
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electron-BP scattering, direct spin-phonon coupling due 
to the phonon-induced strain and the second-order pro- 
cess of electron-BP scattering combined with the hyper- 
fine interaction always has several orders of magnitude 
variation with E at different do and B. However, the 
spin dephasing induced by the hyperfine interaction only 
increases a little with E, which suppresses the large vari- 
ation of the spin dephasing induced by the other three 
mechanisms. Nevertheless, there is still two orders of 
magnitude variation of the spin dephasing when B and do 
are small (Fig Ufa)). The large variation of spin dephas- 
ing induced by the three mechanisms related to electron- 
phonon scattering comes from the fast increase of the en- 
ergy difference A£i2, similar to the analysis of spin relax- 
ation. The spin dephasing induced by the hyperfine inter- 
action is not so sensitive with E. This is because for the 
hyperfine interaction, T% ~ E Z A~ 2 N ', which is obtained 
from the pair-correlation approach, with E z , A and iV 
being the Zeeman splitting energy, the hyperfine interac- 
tion parameter and the nuclear number in the quantum 
dot respectively^ With the increase of the electric field 
E, E z and A keep unchanged and the wave function of 
the ground state is gradually localized on one dot with 
lower potential. This means the effective quantum dot 
size decreases and consequently TV decreases. However, 
as the effective quantum dot size decreases from two dots 
into one dot with the bias field E, N for large E is only 
about a half of that for small E. Therefore, T 2 _1 for large 
field is about two times of that for small field due to the 
decrease of N . However, this increase of T 2 ~ induced by 
the hyperfine interaction is very small compared to that 



induced by other three mechanisms (which have several 
orders of magnitude variation). What should be pointed 
out is that the variation of spin dephasing does not exist 
at T = K. This is because the hyperfine interaction is 
dominant for the spin dephasing at T — K, which keeps 
unchanged with E. 



IV. CONCLUSION 

In conclusion, we propose a scheme to manipulate the 
spin decoherence (both T\ and T<z) in DQD system by 
a small gate voltage. Up to ten orders of magnitude of 
spin relaxation and up to two orders of magnitude of spin 
dephasing can be obtained. To obtain large variation of 
spin decoherence, the inter-dot distance and/or the bar- 
rier hight should be large enough in order to guarantee 
that the two QDs are nearly independent for the small 
bias voltage. At the same time, the effective diameter 
and magnetic field should be small to get as large varia- 
tion as possible. Finally, based on the present available 
experimental technology, the DQD system applied in this 
paper can be realized easily. 11 ' 12 
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